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^ ■ ABSTRACT 

> : 

; Context. A large number of pulsars and X-rays sources are detected in globular clusters. To understand the structure and content of these 
clusters, accurate distances are required. 

Ams. We derive the distances of Terzan 5, Liller 1 and UKS 1 using as a reference a recent distance determination of NGC 6528, based on 
HST/NICMOS and NTT/SOFI infrared photometry. The distance of the metal-poor cluster Terzan 4 was derived from a comparison with M92 
lO in NICMOS bands. 

Methods. Distances of the metal-rich clusters are obtained by comparison of the Horizontal Branch (HB) level of the clusters, relative to the 
reddening line passing through the HB of NGC 6528. We use methods based on NICMOS bands and transformations to J and H magnitudes 
with different assumptions. 

^ Results. Liller 1 and Terzan 4 are found to be at the central bulge distance, UKS 1 is beyond the Galactic center, while Terzan 5 is closer to the 
r Sun than the other four clusters. 

Conclusions. The distance of Terzan 5 is of paramount importance, given the impact of its population of 21 pulsars, which is related to the high 
cluster density. The distance of Terzan 5 is found to be dQ=5.5±0.9 kpc from the Sun, thus closer to us than values given in studies of pulsars 
■ - - in Terzan 5. As a consequence, the higher cluster density is even more favourable for formation of the millisecond pulsars recently detected in 
this cluster. 

Key words. Galaxy: Bulge - Globular Clusters: NGC 6528, Terzan 5, Liller 1, Terzan 4 and UKS 1 - HR diagram 
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1. Introduction 

Globular cluster distances are fundamental to studies of the 
Galaxy structure and stellar population components (e.g. Bica 
et al. 2006). Distances of clusters near the Galactic center are 
particularly interesting, since several of them host a large num- 
ber of X-ray sources and pulsars (Cackett et al. 2006; Ransom 
et al. 2005; Lorimer 2005). Accurate distances are crucial to 
estimate cluster stellar densities, which in turn aflFect stellar in- 
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teraction rates in globular clusters, where dense environments 
favour formation of binaries. Verbunt & Hut (1987; see also 
Hut et al. 1992) reported that 20% of low-mass X-ray binaries 
are found in globular clusters. They also pointed out that Terzan 
5 and Liller 1 might be the densest clusters in the Galaxy. 
Ransom et al. (2005) and Lorimer (2005) reported 21 pulsars in 
Terzan 5. Camilo & Rasio (2005) reported that Terzan 5 and 47 
Tucanae together have 45 millisecond pulsars, as a probable re- 
sult of the "recycling" model (where an accreting neutron star 
is spun up and becomes a fast radio pulsar). This model con- 
nects low mass X-ray binaries (LMXB) to millisecond pulsars 
(MSP). 

The metal-rich globular cluster NGC 6528 is a template ob- 
ject of the bulge population (Bica 1988; Ortolani et al. 1992; 



2 



S.Ortolani et al.: Distances of Terzan 5, Liller 1, UKS 1 and Terzan 4 



Ortolani et al. 1995), and is used here as a reference. The dis- 
tance of NGC 6528 was derived by Momany et al. (2003) from 
a combination of proper motion cleaned optical and infrared 
Colour-Magnitude Diagrams (CMDs), obtained from images 
observed with HST-WFPC2, and the ESQ NTT-SOFI infrared 
camera. The resulting CMDs were compared to fiducial lines 
of 47 Tucanae and NGC 6553 and to Padova isochrones. A 
distance modulus of (m-M)o=14.44+0.06 and a distance from 
the Sun d0=7.7±O.2 kpc were obtained, where the errors were 
derived from a comparison between empirical and a theoreti- 
cal calibration of the HB. The distance errors in the infrared 
are small relative to the optical, because of the minimized red- 
dening dependence. Momany et al.'s cleaned CMDs show un- 
precedented accuracy at the HB level of K(HB) - 13.20+0.05. 
A re-check taking into account calibration uncertainties, indi- 
cates that the distance modulus error should be of the order 
of +0.10. Note that the distance of the Galactic center consid- 
ered previously to be at 8 kpc (Reid 1993) has been recently 
revised to be closer to the Sun at 7.62±0.32 kpc (Eisenhauer 
et al. 2005), 7.5+0.1 kpc (Nishiyama et al. 2006) and 7.2+0.3 
(Bica et al. 2006). 

Ortolani et al. (2001) studied HST-NICMOS Colour- 
Magnitude Diagrams of the reddened bulge globular clusters 
NGC 6528, Terzan 5, Liller 1, Terzan 4 and UKS 1, to derive 
their ages. The results confirmed that these bulge clusters are 
old, with ages comparable with that of 47 Tucanae, as previ- 
ously found by Ortolani et al. (1995). In the present work, we 
concentrate efforts to derive distances for this sample, taking 
NGC 6528 as a reference. 

From the same NICMOS data, Cohn et al. (2002) derived a 
distance of Terzan 5 from the Sun of d0=8.7 kpc. From optical 
NTT-SUSI photometry, Ortolani et al. (1996a) derived do=5.6 
kpc. In the present paper, we readdress the Terzan 5 distance 
issue, given its impact on the cluster density implying in bi- 
nary and pulsar formation. Matsunaga et al. (2005) searched 
for SiO masers in globular clusters and a very luminous one 
was found in the central part of Terzan 5, proposed to result 
from the merging of two AGB stars, likely to occur at this ex- 
treme stellar density. Also, Cohn et al. (2002) reported an RR 
Lyrae in Terzan 5, the first detected at such high metallicities, 
and it was interpreted as due to interaction effects. 

Distances based on optical CMDs were derived for the 
bulge globular clusters Terzan 4 (Ortolani et al. 1997a), Terzan 
5 (Ortolani et al. 1996a), Liller 1 (Ortolani et al. 1996b) and 
UKS 1 (Ortolani et al. 1997b). The purpose of this paper is to 
derive more accurate distances for these clusters making use of 
NICMOS data, as well as to better understand the large num- 
ber of pulsars in Terzan 5. Distances of central clusters are also 
important to establish a locus of survival of old star clusters in 
the central regions of the Galaxy. 

The observations and reductions are briefly reported in 
Sect. 2. In Sect. 3 we discuss cluster metallicities, based on the 
Red Giant Branch (RGB) slope in the infrared. In Sect. 4 cali- 
brations of NICMOS mi 10 and migo photometry are discussed. 
In Sect. 5, distances derived in the present work are presented. 
In Sect. 6 implications for the pulsar rich cluster Terzan 5 are 
developed. Concluding remarks are given in Sect. 7. Details on 
calibrations of NICMOS photometry are given in Appendix A. 
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Fig. 1. Colour-Magnitude Diagrams and mean loci in mno vs. 
mi 10 - mi6o for NGC 6528 and Liller 1 . 

2. Observations and reductions 

Ortolani et al. (2001) observed the sample globular clusters 
with NICMOS on board HST, through the Fl lOW and F160W 
filters, using NICl and NIC2 cameras. The present analysis is 
based mostly on the NIC2 frames. The pixel size of NIC2 is 
0:'075, giving a field of view of 19;'2 x 19;'2 for each frame. 
The description of data, reduction, and a Log of the observa- 
tions are detailed in Ortolani et al. (2001). The typical FWHM 
of the stars was 0!' 1 1 . The zero point PHOTZPT was applied 
to transform instrumental magnitudes into the standard HST 
magnitudes mno and mi6o. 

All the CMDs reach well below the HB. The CMDs of 
NGC 6528 and Terzan 4 reach more than 3 magnitudes be- 
low the turnoff in both the NICl and NIC2-offset pointings. 
The deepest CMD for Terzan 5 reaches about 1 mag. below the 
turnoff". 

3. IVIetallicities of the sample clusters 

Metallicities of globular clusters can be estimated from the 
RGB slope as previously applied to optical and infrared CMDs 
(Ortolani et al. 1990, 1991; Tiede et al. 1997). 

In Figs. [Hand |2] we show the individual CMDs of the four 
metal-rich sample clusters, and corresponding mean loci, plot- 
ted by eye. We estimate an uncertainty of ±0.1 in both mno 
and mno - niieo- In Fig. [3] we also show the metal-poor cluster 
Terzan 4 (Sect. 4.1). 

In Fig. m are plotted the mean loci of the sample clusters, 
by matching their HB level, and subgiant branch (SGB) at the 
HB colour, or the turn-off' in the case of Terzan 4. It appears 
that NGC 6528 is the most metal-rich cluster of the sample, 
closely followed by Terzan 5, and next by Liller 1 and UKS 1, 
in this order. The RGB of Terzan 4 is steeper due to its low 
metallicity. The RGB slope measured in the CMD of NGC 
6528 is Amiio/A(miio - mieo) = 11.0, while for UKS 1 it is 
Ami io/A(mi lo -mi6o)= 15.0+1 .0, the latter value being the most 
uncertain one. 
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Fig. 2. Colour-Magnitude Diagrams and mean loci in mi lo vs. 
mno - mi6o for Terzan 5 and UKS 1. 
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Fig. 3. Colour-Magnitude Diagram of Terzan 4, and mean loci 
of M92 overplotted, in mno vs. mno - nii6o- 



Recent high resolution spectroscopic determinations con- 
firm this ranking for NGC 6528 (Zoccali et al. 2004), Terzan 
5 and Liller 1 (Origlia et al. 2002; Origha & Rich 2004), and 
UKS 1 (Origha et al. 2005). 

Reddening and metallicity values for the sam- 
ple clusters revised by Barbuy et al. (1998) and 
the compilation by Harris (1996), as given at 
'http://www.physics.mcmaster.ca/Globular.html, as well as 
recent spectroscopic determinations are shown in Table 1 . 

4. Distance from NICMOS instrumental 
magnitudes 

For the derivation of distance, two methods are possible, one 
using comparisons with template clusters in the instrumental 
NICMOS bands, and the second by transforming them into J 
and H. 

Schlegel et al. (1998)'s reddening law (1998) was em- 
ployed to obtain the infrared extinction slope Aiio/(Aiio-Ai6o) 
in the HST/NICMOS bands: Aii()/E(B-V)=1.186, Ai6o/E(B- 
V)=0.634, and Aiio/(Aiio-Ai6()) = 2.148. As a check, we also 
calculated the slope using the Lee et al. (2001) NICMOS ex- 
tinction law based on synthetic spectra, obtaining a similar 
value Aiio/(miio-mi6()) = 2.03. 

Fig. |5] shows the observed mean loci for the 5 clusters, in 
NICMOS mno vs. mno-miso- The reddening line Aiio/E(miio- 
mi6o) = 2.148 passing through the HB of NGC 6528 is also 
shown in Fig.|5] Were all the clusters at the same distance, then 
the HB of all of them would be on the same line. Therefore, 
the mno difference between each cluster HB and the reddening 
line represents the distance modulus difference between NGC 
6528 and the cluster, assuming that they share the same metal- 
licity (or same intrinsic SGB colour and HB absolute value). 
The advantage of this method is that reddening from NICMOS 
colours is not needed. 

Fig.|5]shows that Liller 1 and NGC 6528 differ significantly 
in reddening. Terzan 5 is located above the reddening line, thus 



at a shorter distance, whereas UKS 1 is below the line, and 
therefore farther than NGC 6528. The distances derived are 
given in column 5 of Table 2. 

4.1. Distance of the metal-poor cluster Terzan 4 

Given the difference in Horizontal Branch morphology be- 
tween metal-poor and metal-rich clusters, the method described 
in Sect. 4 cannot be applied to Terzan 4. The distance of 
Terzan 4 was derived by comparing its CMD to the mean lo- 
cus of M92 (Lee et al. 2001) observed in the same NICMOS 
bands (Fig. [3j. Lee et al.'s reddening law in the NICMOS 
bands was adopted. The best fit gives A(miio-mi6o)(Tz4-M92) 
= 0.73+0.02 and A(miio)(Tz4-M92) = 1.6+0.1. Adopting E(B- 
V) = 1.87 E(miio-mi6o) = 1.38 and Aiio/E(B-V)=1.16 (Lee et 
al. 2001), we get Ano = 1.38 x 1.16 = 1.60 and A(m-M)o(Tz4- 
M92) = 0.0. Finally Ay = 1.38 x 3.3 = 4.5 mag. 

This means that Terzan 4 is at the same distance as M 92, 
which is 8 kpc according to Harris (1996). Therefore the dis- 
tance of Terzan 4 is R0=8.O±O.3 kpc. This is very close to 
the Galactic center Its extinction is rather below the 7 magni- 
tudes we derived in the optical, but the distance is compatible 
with our previous estimations of 8.3+0.7 kpc (Ortolani et al. 
1997a) and 7.3 kpc (Barbuy et al. 1998). The present determi- 
nation gives the most reliable distance and reddening values for 
Terzan 4. 

4.2. Errors 

Distance errors in the method described above are dominated 
by uncertainties on: 

(a) The measurements of SGB colours and luminosities at 
the level of the HB. Differential reddening across the images 
is the main source of errors, and such effects are common to 
all calibrations. Differential reddening together with some field 
contamination of the lower Galactic latitude clusters (such as 
Liller 1) are the main causes of scattering in the CMDs. The 
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intrinsic photometric errors are negligible because of the very 
high signal to noise ratio at the HB level, given that the pho- 
tometry reaches below the turnoff. In the CMDs they are of the 
order of 0.07 in colour and 0.15 in magnitude. These values 
were obtained by measuring the HB in different parts of the 
images and/or comparing NICl and NIC2 cameras. However, 
since the CMD scatter is mainly in the direction of the red- 
dening line, the effect propagated on the distance is reduced to 
about 0. 1 mag. A metallicity uncertainty of about 0.2dex would 
imply the same error (0.1 mag). 

(b) The slope of the reddening law. The formal error on 
the distances introduced by the uncertainty on the slope of the 
reddening line (Schlegel et al. 1998; Fitzpatrick 1999) is con- 
siderably lower, of the order of 0.02 mag. in the full colour 
range. Second order effects may arise due to the wide passband 
of the FUOW and F160W filters (the first is about twice that 
of the standard J, and the second 35% wider than the H band). 
Possible non linear effects due to the distorted passbands may 
be an important uncertainty of this calibration and cannot be 
easily quantified. While the second order reddening effects are 
well known in the optical photometry, in the infrared the effects 
on the standard K band have been only recently inferred (Kim 
et al. 2005). 

5. Reddening and distance from NICI\/IOS 
calibrated data 

J and H calibrated magnitudes and colours have been obtained 
by adopting linear transformation equations by Stephens et al. 
(2000) and Cohn et al. (2002). Both authors used the NICMOS 
Data Handbook standard stars, whereas Stephens et al. added 
Baade's Window giants observed both with NICMOS and from 
ground in JHK. More details on the calibrations are given in the 
Appendix. 

In Fig. |6]we show the J level of the HB ws. J - H of the 
SGB at the HB level (SGB@HB) for the 4 metal-rich sample 
clusters. In Figs.|5]and|6] the distance determination is indepen- 
dent of reddening values and zero point calibration. A source of 
uncertainty is the choice of the reddening law, but differences 
between them in the infrared are small (Schlegel et al. 1998; 
Fitzpatrick 1999). Moreover, there is a general agreement in 
the literature that, in the infrared, the reddening slope is con- 
stant and independent of the dust characteristics. 

Transformations of E(J-H) into E(B-V) are reported in 
the Appendix. Assuming for NGC 6528 the reddening E(B- 
V)=0.46 (Zoccali et al. 2004), we get for Terzan 5, Liller 1 and 
UKS 1 the AE(B-V) values. In order to obtain the total redden- 
ing, that of NGC 6528 itself has to be added, and the results 
are given in Table 2 (columns 9 and 13 for the two calibrations 
respectively). A good agreement is found with optical CMDs 
(Table 1) for Terzan 5. For Liller 1 and UKS 1, Cohn et al.'s 
calibration is more compatible with the literature. 

In Table 2 are given distance determinations for the sam- 
ple clusters, assuming for NGC 6528 the value d^-l .1 kpc 
from Momany et al. (2003), and using the independent cali- 
brations from instrumental to J and H bands by Stephens et al. 
(2000) and Cohn et al. (2002). Columns 8 and 12 of Table 2 
give the distance modulus difference between each cluster and 



NGC 6528, measured along the reddening line (Figs. |5] and |6]l, 
for each calibration. This makes this procedure reddening inde- 
pendent. The differences in distances between the two calibra- 
tions arise from the different calibration coefficients, and not 
from photometric zero points. 

In Table 2 we calculated the distances of UKS 1 and Liller 
1 for two metallicity values: the same as that of NGC 6528 
([Fe/H]=-0. 1), and a correction of A[Fe/H]=0.5, making use of 
Padova isochrones (Girardi et al. 2000). The colour changes by 
A(J-H)«0.1, corresponding to a decrease of 0.25 magnitudes in 
the distance modulus. 

Terzan 5 distance estimates of 4.6 and 5.6 kpc, put it closer 
than the Galactic center, whereas a mean of 6.5 and 9.1 kpc 
for Liller 1 place the cluster near the Galactic center UKS 1 
estimates of 9.3 to 12.9 kpc suggest that the cluster is beyond 
the Galactic center. 

6. Pulsars and density of Terzan 5 

Terzan 5 distance estimates range from 4.6 to 6.3 kpc (Table 2), 
shorter than Cohn et al. (2002)'s value of 8.7 kpc. Cohn et al. 
derived this distance in two ways: (i) One using a comparison 
of HB level of NGC 6528 in the H band, adopting H(HB)=13.1 
from Davidge (2000). However it can be checked in Davidge's 
Fig. 5 that instead H(HB)=13.5 for NGC 6528, and the same 
value is found in Momany et al. (2003). Terzan 5 H(HB) is also 
different: Cohn et al. adopt H(HB)=13.0 (in fact incompatible 
with his Fig. 6 that indicates H(HB)=14.0) whereas we have 
H(HB)=13.76 and 13.9 according to Stephens et al. or Cohn 
et al.'s calibrations, (ii) Cohn et al. used the Baade's Window 
Red Giant Branch as reference, and deduced E(B-V)=2.18 for 
Terzan 5. Using the distance from Ortolani et al. (2001) of 5.6 
kpc, together with their lower reddening, a distance of 8.7 kpc 
was found for Terzan 5 by Cohn et al. 

In conclusion, the incorrect H(HB) values in (i), and a low 
reddening of (E(B-V)=2.18 for Terzan 5 in (ii), may explain 
Cohn et al.'s different distance estimates relative to ours. 

In the present work (Table 2), estimated distances of Terzan 
5 are 6.3, 4.6 and 5.6 kpc. The three values are consistent, rein- 
forcing the results, given that the first method is independent of 
photometric calibrations. An average distance d© = 5.5+0.9kpc 
is adopted for Terzan 5. 

Camilo & Rasio (2005) presented evidence that globular 
clusters with pulsars have densities p > 1000 Lo/pc^. In their 
Fig. 6 Terzan 5 is reported to be one of the densest clusters, 
with p > 10^'^ Lo/pc^, assuming however a distance of 7.6 kpc 
from the Sun (Harris 1996). 

A closer distance causes a higher calculated stellar den- 
sity (Lo/pc-^). Luminosity scales with R^, whereas the volume 
scales with R^, therefore the density increases for shorter dis- 
tances. The present results imply that the density would in- 
crease by 30-40%. 

7. Conclusions 

Using NICMOS on board the Hubble Space Telescope 
(Ortolani et al. 2001), we derive relative distances for Terzan 5, 
UKS 1 and Liller 1, based on the distance value of NGC 6528 
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Table 1. Literature reddening and metallicity (Barbuy et al. 1998; Harris 1996 and high resolution spectroscopy by 1 Origlia & 
Rich 2004; 2 Carretta et al. 2001; 3 ZoccaU et al. 2004; 4 Origlia et al. 2005; 5 Origlia et al. 2002) 









Barbuy et al. 1998 


Harris' compilation 
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1 


b 


E(B-V) 


[Fe/H] 


E(B-V) [Fe/H] 


[Fe/H] spec 


ref. 


NGC 6528 


1°14 


-4° 18 


0.52 


-0.2 


0.56 -0.17 


-0.1 


2,3 


Terzan 5 


3°81 


1°67 


2.39 


0.0 


2.37 -0.28 


-0.21 


1 


Liller 1 


354°81 


-0°16 


3.0 


+0.2 


3.00 +0.22 


-0.3 


5 


UKS 1 


5°13 


0°76 


3.1 


-0.5 


3.09 -0.5 


-0.78 


4 


Terzan 4 


356°02 


1°31 


2.31 


-2.0: 


2.35 -1.60 


-1.6 


1 



Table 2. Horizontal branch magnitudes, colours and distances for different photometric calibrations. The mno and J magnitudes 
refer to the HB, colours refer to the SGB at the HB level. Amno and AJ are the distance moduli differences relative to the 
reference cluster NGC 6528. Assuming for UKS 1 and Liller 1 a metallicity of" the same as NGC 6528 of [Fe/H] =-0.1; * more 
probable value of [Fe/H] =-0.6. 



NICMOS instrumental Stephens et al. calibration Cohn et al. calibration 
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d(kpc) 


J(HB) 


J-H(HB) 


AJ 


E(B-V) 


d(kpc) 


J(HB) 


AJ 


E(B-V) 


d(kpc) 


NGC 6528 


-0.36 


16.71 




7.7 


14.2 


0.58 




0.46 


7.7 


14.11 




0.46 


7.7 


Terzan 5 


0.44 


17.99 


-0.44 


6.3 


15.22 


1.19 


-1.10 


2.82 


4.6 


15.09 


-0.70 


2.33 


5.6 


Liller 1° 


0.75 


16.66 


0.57 


10.0 


16.79 


1.42 


-0.36 


3.74 


6.5 


16.79 


0.36 


3.04 


9.1 


Liller 1* 








8.9 










5.8 








8.1 


UKS 1° 


0.70 


20.56 


1.57 


15.9 


17.69 


1.39 


0.68 


3.59 


10.6 


17.69 


1.34 


2.94 


14.3 


UKS 1* 








14.2 










9.3 








12.9 



derived by Momany et al. (2003). The distance of the metal- 
poor cluster Terzan 4 was derived from a comparison with M92 
in the same NICMOS passbands. 

From the diagram /(HB) vs. J - //(SGB@HB) (Fig. ill, 
Terzan 5 is found to be closer to the Sun than NGC 6528, 
whereas Liller 1 is at a comparable distance, and UKS 1 ap- 
pears to be more distant. Assuming a distance to the Galactic 
center of Rgc=8 kpc (Reid 1993), or closer distances of 
Rgc=7.2 to 7.6 kpc (Eisenhauer et al. 2005; Nishiyama et 
al. 2006; Bica et al. 2006), NGC 6528 is essentially at the 
Galactic center distance. Liller 1 could be somewhat beyond 
the Galactic center 

The relative metallicity of the sample metal-rich clusters, 
based on RGB slope in the infrared, show that NGC 6528 is 
the most metal-rich cluster, closely followed by Terzan 5, Liller 
1 and UKS 1 . Terzan 4 is clearly more metal-poor. The lower 
metallicity of UKS 1 favours the shorter distance of 9.3 kpc 
based on Stephens et al. (2000)'s calibration, placing it in the 
bulge volume. 

We have shown that Terzan 5 is at d0=5.5 kpc from the Sun, 
thus significantly closer than the value of do=8.7 kpc calculated 
by Cohn et al. (2002), or d0=7.6 kpc adopted by Camilo & 
Rasio (2005). Recently, Ransom (2006) pointed out that if the 
distance of Terzan 5 were really 8.7 kpc, then the current mod- 
els of electron distribution in the Galaxy would overestimate 
the integrated electron density towards the cluster. For 8.7 kpc, 



models predict a dispersion measure of 530 pc/cm^, while for 
the pulsars in Terzan 5, the dispersion measured is estimated to 
be of 239 pc/cm^. This discrepancy is better explained by the 
distance of 5.5 kpc. The present shorter distance has also im- 
portant implications concerning the cluster density. A 30-40% 
decrease in distance using Camilo & Rasio's (2005) value leads 
to a density increase by the same fraction, which is important 
for the calculation of the stellar collision rate, leading to the 
formation of low mass X-ray binaries and millisecond pulsars 
in this cluster (Edmonds et al. 2001; Ransom et al. 2005). 

Appendix A: JH calibrations of NICIVIOS 
photometry 

The calibrations we analysed in detail are: 

a) Stephens et al.(2000): 

J = -0.344(miio - mi6o) + const 
J - H = 0.96(miio - mi5()) + const 

b) Cohn etal. (2002): 

J = -0.335(mno - mieo) + const 
J - H = 0.68(mno - mif,o) + const 

c) Girardi etal. (2000) for Z^0.02: 

J = -0.32(mno - mieo) + const 
J - H = 0.77(mno - mif,o) + const 

d) Girardi et al. (2000) for Z^O.0004: 

J = -0.2(miio - mi6o) + const 
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Fig. 4. Mean loci in mno vs. mno - migo for the five sam- 
ple clusters, by matching relative to NGC 6528 their HBs and 
SGB@HB, or the turn-off in the case of Terzan 4. NGC 6528: 
black circles Terzan 5: open triangles; UKS 1: black triangles; 
Liller 1: open pentagon; Terzan 4: black square. 




-0.5 0.5 1 



Fig. 5. Mean loci in mno vs. mno - mi6o for the five sam- 
ple clusters, in ST magnitudes. The solid fine is the reddening 
line with a slope Aiio/E(mno-i6o)=2.15, passing through NGC 
6528. NGC 6528: open circles; Terzan 4: open squares; Terzan 
5: open triangles; Liller 1: open circles; UKS 1: open stars. 



Fig. 6. J level of the HE vs. J - H of the SGB at the HB 
level for the four metal-rich clusters, derived with two differ- 
ent transformation equations. The solid line is the reddening 
Une with a slope Aj/E(J-H)=2.76, passing through NGC 6528. 
Symbols: open triangles: transformation equations by Stephens 
et al. (2000); filled triangles: transformation equations by Cohn 
et al. (2002). 

J - H = 0.86(miio - mi6o) + const 

A.1. Uncertainties in tlie colour transformation 

The errors on the distances and reddening values are domi- 
nated by the transformation coefficients (colour terms), from 
the NICMOS instrumental system, into the JH magnitudes. 
More reddened clusters are more dependent on the redden- 
ing/caUbration slopes, which is the case of UKS 1 and Liller 
1. 

The calibration of wide band photometry in non-standard 
systems, such as NICMOS, is not trivial. The estimation of un- 
certainties has no simple solution because standard stars are in 
general not available for extreme temperature and reddening 
values. Non-linearities and second order effects may be non 
negligible. In the following we analyse separately the two dif- 
ferent calibrations that we used to derive the distances (Table 
2), and compare them with other caUbrations in the literature. 

A.2. Stephens et ai.'s calibration 

Stephens et al.'s calibration (equations a) has a J slope coeffi- 
cient (-0.344) similar to Cohn et al. (2002)'s calibrations, but 
their H slope is larger (-0.305, vs. ^ -0.1). This provides a 
flat slope for the J-H transformation equation, and somewhat 
redder colours for very red stars. This explains the relatively 
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Fig. A.l. Colour calibrations from Stephens et al. (2000) (solid 
line), Cohn et al. (2002) (dotted) and Origlia & Leitherer 
(2000) (dashed). The mno-mieo colours of NGC 6528, Terzan 
5 and Liller 1 are shown. 

shorter distance scale for the most reddened clusters (Liller 1 , 
UKS 1), combined with high reddening. 

Stephens et al. (2000) checked fits with simulations using 
infrared spectra of reddened stars. Assuming that bulge field 
giants are similar to cluster giants, this method potentially of- 
fers the most reliable transformations for bulge clusters. On the 
other hand, data for Baade's Window stars show some scatter, 
and they cover a limited colour range, producing errors of ±0.1 
on the colour slope and ±0.08 mag. on the zero point, and an er- 
ror up to ±0.2 mag at the extremes of the colours of the sample 
clusters. 

A.3. Cohn et al.'s calibration 

The calibration used by Cohn et al. (2002) (equations b) is 
based on the NICMOS five solar neighbourhood standard stars, 
distributed over a wide temperature/colour range. The J-H 
colour term is lower than in Stephens et al.'s calibration (0.76 
vs. 0.96), resulting in bluer colours and increasing distances for 
more reddened clusters (Table 2). The formal errors, from the 
linear interpolation, are similar to those in Stephens et al. They 
are marginally compatible with the high metallicity of Liller 1 
and UKS 1, and with Liller 1 JHK photometry (Frogel et al. 
1995). The JH CMDs of NGC 6528 calibrated with Stephens 
et al. and Cohn et al.'s equations have been compared with 
Momany et al. (2003)'s ground based photometry. The colours 
were found in better agreement with Stephens et al.'s calibra- 
tion, whereas Cohn et al.'s calibration gave somewhat bluer 
colours. 



A.4. Other calibrations 

Cohn et al.'s calibration equations are similar to others in the 
Uterature, such as Macri et al.'s (2001). In the latter case the 
calibration was derived from solar metallicity Kurucz models, 
convolved with NICMOS passbands and compared to the stan- 
dard JH system. Reddening effects were not considered. They 
derived small formal errors, one order of magnitude less than 
those of Cohn et al. (2002). 

A colour calibration was obtained by Origlia & Leitherer 
(2000). They used Bessell et al. (1988) model atmospheres 
to transform NICMOS colours into the JHK system, and con- 
cluded that the transformations are gravity and metalUcity in- 
dependent, but a cubic polynomial relation is required to fit the 
data. The efi'ect of such an interpolation is that for NGC 6528 
and Terzan 5, with reddening values lower than UKS 1 and 
Liller 1, the results are similar to Cohn et al. (2002) and Macri 
et al. (2001) cahbrations, but the colours are considerably bluer 
for more reddened clusters, implying very large distances for 
UKS 1 and Liller 1 . The effect of linearity deviations is seen 
in their Fig. 1 1, as well as a discrepancy with the reddest stan- 
dard stars. Fig. A.l shows the different colour calibrations in 
the metallicity range of our clusters. 

From the analyses above we conclude that systematic errors 
of the order of ±0.2-0.3 mag., related mostly to colour transfor- 
mation slopes, affect distance moduU, in particular those of the 
most reddened clusters. 

A.5. Calibrations with Padova stellar evolution models 

In order to clarify the nature of the discrepancies between the 
described calibrations, we performed a further independent test 
based on the Girardi et al. (2002) NICMOS isochrones in the 
mi 10 and miso bands. The authors derived the isochrones from 
ATLAS9 models (CasteUi et al. 1997) and a NICMOS trans- 
mission system provided by the NICMOS staff. Comparing 
these isochrones with those by Girardi et al. (2000) in the JHK 
system, we obtained the transformation equations c) and d) 
for metallicities of Z=0.019 (solar), and a low metallicity of 
Z=0.0004, in the full colour range of the isochrones (about 
0.2<miio-mi6o<l.l). For solar metallicity the linear interpo- 
lation gives a colour term of -0.32, somewhat lower than that 
obtained by Stephens et al. (2000). At lower metallicities the 
coefficient decreases to -0.2 at Z=0.0004. The colour coeffi- 
cients in this calibration are intermediate between Stephens et 
al. and Cohn et al.'s. The result of this test indicates that the 
metallicity is not the major source of the discrepancy between 
Stephens et al. (2000) and the synthetic photometry. We ar- 
gue that an important source of discrepancy could be redden- 
ing, in agreement with Lee et al. (2001). Using their redden- 
ing law calculated for NICMOS bands they derived an average 
A(miio)/E(B-V)=1.16 (for SGB and RGB stars), while Rieke 
and Lebofsky (1985) give Ay/E(B-V)=0.864. A star reddened 
by one magnitude in the mno-migo colour will imply A(miio)- 
Aj=0.55 mag. It is evident that the colour terms coming from 
synthetic photometry are by far too small to account for the 
reddening effects. 
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We conclude that a correct calibration cannot be obtained if 
the reddening is not known in advance. A solution could be an 
iterative process deriving a first guess reddening and proceed- 
ing further with subtraction of the reddening and calibration of 
the dereddened data, adding the obtained final reddening to the 
calibrated data. 
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